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Abstract 


A charging n odfel is dfevelbped for geosynchronbus, three-axis stabilized 
spacecraft when i nder the inHuence of a geomagnetic sUbstorm. The differential 
charging potentia .s between the thermally coated or blanketed outer surfaces and 
metallic struCtun of a spacecraft are determined when the spacecraft is immersed 
in a dense plasma cloud of energetic particles. The spacecraft -to-environment 
interaction is determined by representing the charged particle environment by 
equivalent current source forcing functions and by representing the spacecraft by 
its electrically equivalent circuit with respect to the plasma charging phenomenon. 
The charging mbdel includes a sun/earth/spacecraft orbit model that simulates the 
sun illumination conditions of the spacecraft outer surfaces throughout the orbital 
flight on a diurnal as well as a seasonal basis. Transient and steady-state numer- 
ical results for a three-axis stabilized spacecraft are preserlted. 


This work «as performed under a General Electric !5pace Division 197G Internal 
Research and Development Program, No. 768084161. 


I. iNrftmMT.TiON 

from thO geosynehratious orbitihg NASA satellitfes ATS-5 and 
ATS-e • • has indicated that the surface of these satellites can charge to hundreds 
of volts when in sunlight and thousands of volts (up to -10 kV) when in eclipse. 

Data transmitted from these satellites during these charging events has indicated 
the existence of transient fluxes of energetic particles. It has been suggested^ 
that these clouds of energetic particles are injected into the IjcaLmidnight-to- 
dawn region of the geosynchronous altitude during geomagnetic substorm activity. 
Consequently, during a geomagnetic substorm, spacecraft at altitudes greater than- 
three Earth radii, in the- local time sector from just before midnight to past dawn, 
occasionally will encounter and be immersed in.a dense plasma cloud of energetic 
particles. . it has been further postulated that this charged particle environment is 
the major cause of spacecraft charging. That is, in the steady-state, every isolated 
part of a spacecraft immersed in the space environmental plasma will come into 
electrical 6<iuilibrium by developing surface charges of the proper sign and magni- 
tude such that the net current - represented by the deposition and release of 
charged particles from the surface of the spacecraft - is zero. The equilibrium 
potential of the surface of the spacecraft is the potential difference between the 
surface and ambient plasma sheath. The most important contributors to the equil- 
ibration currents are the primary plasma electron artd proton arrivals at the sur- 
face and the photoelectrons released when sunlight illuminates the surface, in 
addition, the Contributions of secondary electrons released from the surface Under 
primary proton or electron impact and possible electron reattraction to the surface, 
are also significant and must be considered in a complete analysis of the problem. 

In the coincidence with the geomagnetic substorm activity in this local time 
quadrant, is the occurrence of anomalous events on-board satellites in geosyn- 
chronous orbit when immersed in the substorm plasma. Specifically, anomalous 
behavior experienced by several satellites has included^ control circuit switching 
power system failure, sensor data noise, thermal Control degradation, and tele- ’ 
metry logic switching. There is a growing body of evidence which demonstrates 
the dependence of satellite anomalous behavior on geophysical parameters such as 
local time and geomagnetic activity. 

Consequently, it has been postulated*' ^ that the anomalous behavior of syn- 
chronous spacecraft is due to electrostatic charging of the various spacecraft 
surfaces to large negative potentials and their subsequent discharging. The elec 
tromagnetic pulses produced by the discharges contain enough energy to interact 
with electronic logic circuits at distances of tens of centimeters, and cause Voltage 
spikes large enough to change logic steles, other data from spacecraft indicate 
that repeated discharging also results in the degradation of thermal control surfaces 
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1 llu>, II la Ihe differential charelnB of the vaMous Itiermally blankelod or coaled 
o or aorfaco n.a.cHaU nllh reapco, ,o each olhcr and .llh reapoo. lo Ihe apace- 

ro„ Ta a ' when Ihe dleleclrlc 

Llano! t '* that la one of the major cauaos of 

satdliit^ aiiornaloUs behavior, 

which adreli'Tle '*"** **^*^'' ‘develop a apacecraft charging almulatlon model 

snacecran ^■’““'f.t.-.o-cnvlronmen. Inlerac.lon when Ih. 

apacecraft immeraed In .he charged particle cnvironmenl that 1. encunlered a, 

g oaynchronoua orbit during a geomagnetic aubstorm. Purther, the principal 

polan lata between Ihe thermally blanketed or coated outer aurfaeea and the metallic 
atreelure of the apacecraft. To determine the apacecraft -to-envlronmenl Inter- 
action, the equivalent electrical model of the apacecr.n with reaped lo the charg- 
ing phenomenon la developed and the plaama environment la repreaentod by equlva- 

oartZ 'h^ charged 

and th! aurrerrr "" 

current u a«ondary electron emlaalona aa dlaeharglng 

nt sources. The spacecraft outer surface configuration is represented by 

• , ^ geometrical model and a solar/eai^th orbital model are 

alao developed lo determine the aun-.llumtnatlon condition of the outer aurfaeea aa 
^function of ap.cec.1, orbital poaltlon. That la, the geometrical and orbital 
models are used to determine whether a surface is sun -illuminated self -had d 
or earth-ahadowed. in ad-lltlon. when a aurface la aun-lllumlZd th mtenar ’ 

of the lllumina ton, which la a function of the aun/apacecratt aurface aapect aZ 
IS also determined by the models, 

mode?^’'®'? ’ However the 

niodela are baaed primarily on the mialyala of apln-alabillred apacecraft. Fo^ Ihe 

particular types of spacecraft analyacd, there was conalderable aeaaonal and 
lurnal variation of the exposed metallic are. Illuminated by the aim. However 
due o the restricted location of thermal blanket malarlala and Ihe exlcrnal atruc- 
mra, orm of spln-slabltlxed sp.cecr.fi. ihere were small diurnal varlatloi! in 
h the amolml and location of the thermal blanket material areas tllummated bv 
the sun during the mldhlghl-to-dawn local lime quadrant r„n 
was made to determine the sun-lllumlnatlon conltoZ'exZTdZtiraZ^^^^^^ 

dtuZl'!,' ^elusZ' 1 'he 

dleleeirL and "’7''“'’'”" of 'he amount and location of holh Ihe 

The moW dl d! T‘"' """""■e'e'* »y the sun are considerable, 

condlllon of alUMh ” ' eelcrmlnes the variations of sun-lllumlnallon 

of the exposed surfaees throughout the dally orbital path in addition 
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to considtsrlne the seasonal changes. Also, the previous models have hot con- 
sidered the intrinsic capacitance of the spacecraft structure vHth respect to the 
plasma sheath. In this paper, the structural capacitance has been included In the 
simulation model, and it will be shown that the structural capacitance has a signif- 
icant influence on the transient response. 

In the following sections, the development of the spacecraft charging model is 
discussed in general terms and includes a discussion of the plasma model, the 
electrical model, the geometrical model, a solar/earth orbital model, the mate- 
rial properties and configuration definition, and the numerical integration 
approach. Spacecraft charging results are given for a geosynchronous satellite 
during the midnight -to -dawh local time quadrant for the fall-equinox artd winter- 
solstice seasonal periods. 


1 Sl‘ ACKCIIAFT CllAHOINt; MOIIKL DFN KI.Ol’MFNT 

The spacecraft charging model development can best be described in terms of 
the flow chart shown in Figure 1. The S/C charging model consists of four sepa- 
rate models: a plasma model, an electrical model, a S/C geometrical model; and 
a solar /earth orbital model. The plasma model represents the charging and dis- 
charging mechanism of the ambient plasma with respect to the spacecraft by equiv- 
alent current sources. The current sources, which are dependent on the particle 
energy distribution functions, constitute the forcing functions of the charging model 
equations. The electrical model defines the lumped element equivalent circuit 
representation of the spacecraft surfaces with respect to the electrostatic charg- 
ing phenomenon. The plasma model and electrical model are combined to form 
the nonlinear spacecraft charging equations. The spacecraft geometrical model 
defines the spacecraft outer surfaces in terms of approximate planar surfaces and 
curved surface projections and defines the vertices of all planar and curved 
surfaces in terms of a spacecraft reference coordinate system. The Solar/earlh 
orbital model determines the location of the spacecraft with respect to the sun and 
the earth. The geometrical model and the solar/earth orbital model are combined 
to determine the variation of the sun -illumination conditions of the outer surfaces 
with respect to orbital position. 

To complete the modeling, the surface material properties and configuration 
are defined. The surface material properties that are most important in a space- 
craft charging analysis are: the relative dielectric constant, the variation of the 
surface resistivity with respect to clecirical stress level, and the variation of the 
bulk resistivity with respect: to electrical stress level. The material configuration 
definition describes the location of the various thermal blanket and surface coating 
materials. 



Figure 1. Spacecraft Charging Model Flow ChaM 


The elements of the flol» chirt ivill now be discussed in greater detail. 

2. 1 Plasifiii Model 

As discussed previously, spacecraft at geosynchronous orbit will occasionally 
eneountet- energetic charged particle fluxes artd these fluxes will cause the various 
outer surfaces to charge to large potentials. The charging and discharging 
mechanism of the ambient plasma with respect to the spacecraft cart be represented 
of simulated by equivalent curreht sources that become the forcing functions of the 
charging model equations. The constituent particle fluxes that affect the charging 
of a surface afe pfotorts, electrons, photoelectrons, and secondary electrons. 

Several investigators®' ®' ® have approximated the energy distribution of the 
pdrtlcle fluxes as measured on ATS-5 by a Maxwell-Soltzmanrt (M-Bl energy dis- 
tribution and, further , have assumed the particle fluxes to have an omnidirectional 
energy distribution; that is, the energy distribution of the particles is identical in 




each dlfdctloh. thus, it Is possible to chafactcrl^.e the partible energy distribu- 
tions by a thermal energy (Kr) expressed In electron volts. This approximation is 
used in this paper to simplify the analysis and to provide an insight into the problem 
that might otherwise be obscured by a more complex approach. 

The total current flowing into the outer surface of a spacecraft is 


where A is the surface area and J.j. is the total positive current density into the 
surface and is given by 

Where J is the incident proton current density, is the secondary electron 
current density produced by incident protons, is the incident electron current 
density, is the secondary electron current density produced by incident elec- 
trons, and Jpjj is the photoelectron current density. 

A charged surface at a given potential in a charged particle environment will 
accelerate particles cf the opposite polarity artd repel particles of the same 
polarity. Thus, assuming an omnidirectional Wlaxwelllan energy alstributloh, the^ 
fraction of ambient plasma electrons reaching a large surface at a potential V is 


exp 

® ^o 



V<0 


(3) 


where N is the incident electron density, N is the ambient electron density, V 

© “q 

is the potential of the surface under consideration, K is Boltzmann's constant, e 
is the charge of an electron, and is the absolute temperature of the M-B elec- 
tron energy distribution. 

The average ambient electron current density incident to a neutral surface is 
given by 


where is the average ambient electron current density and v^ is the mean 
ambient Aermal velocity. 

Thus, from Eqs. (3) and ( 4 ), the average electron current density incident to 
a large surface at potential V is 
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Je = Jfi ^ • V <0 

® ®o 


where 


T = — - 
e e 


( 3 ) 


( 6 ) 
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id the equivalent temperature, expressed in Volts, of the distribution approx- 
imating the plasma eleetron energy distribution, and is the eleetron current 
density incident to a large surface at potential V. For a positive surface 

= , V>0 (g) 

that is, a surface at a positive potential will attra-!t oppositely charged particles 
but cannot extract more particles from the plasma environment than the ambient 
particle density Similarly, the proton current density incident to a large 

surface at potential V is 

Jp-Jp^e P . V>0 (9) 

and 

jp = Jp^ , VsiJ (10) 

where ^Po is the average ambient proton current density incident to a neutral sur- 
face, and Tp is the equivalent temperature of the M-6 distribution approximating 
the plasma protoh energy distribution add is expressed in volts. 

In addition to the above charged particle fluxes, there will be secondary emis- 
sion electrons as well as photoelectroh emissions. !&oth types of charged particles 
will be repelled by a surface at a negative potential and attracted by a surface at a 
positive potential. (!!onsequently, based on the previous discussion, the secondary 
electrons leaving a surface of potential V is given by 
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( 12 ) 


-V/T^ 

Jsp=Jsp « " . V^O 


Jsp = Jsp^ . VsO 


( 13 ) 


( 14 ) 


whfefe and are the average secondary electron current densities emitted 
^rom a neutral surface produced by incident electrons and protons, respectively. 
Jg^ and Jsp are the secondary electron current densities emitted by the incident* 
electrons and protons, respectively, and is the eduivalent temperature of the 
M-B distribution representing the energy distribution of the secondary emission 
electrons and is expressed in volts. The secondary emission electron current 
densities are directly related to the Incident particle current densities. It is 
assumed that the secondary electrons emitted from a neutral sUfface are related 
to the incident pafticles by a fixed Constant and can be expressed as 



Je^e 


( 15 ) 



= L 

P P 


( 16 ) 


Inhere fg is the ratio of secondary electrons to incident electrons and f is the ratio 
of secondary electrons to Incident protons. In general, the secondary emission 
constants, f^ and fp, arill have different values for dielectric and metallic surfaces. 

A similar development holds for the photoelectron emissions produced by sun 
Illumination. The photoelectron current is directly proportional to the Intensity 
of sunlight which is related to the angle of incidence. Consetjuently, the photo- 
electron current density emitted from a surface at potential V can be expressed as 


^ph ■ "^ph 


■V/T, 


ph 


cos a 


V >0 


( 17 ) 
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"^ph 



cos a 


V sO 


<1H) 


where J < Is the photoelectron current density emitted from an illuminated surface 
at potential V, ts the average photoclectron current density emitted from an 

iliuminated neutral surface, Tpj^ is the equivalent temperature of the M-H distri- 
bution representing the energy distribution of the photoelcctrons expressed in volts, 
a Is the angle between the sun-line and the surface normal vector (sun/spacccraft 
surface aspect angle) and 


cos a 


cos a for lo| < jr/2 

0 for Id s ?r/2 (self-shadowing conditions) 


(ly) 


The total positive current density into a surface can take one of four possible 
forma depending on the polarity, positive or negative, of the surface potential and 
the presence or absence of sun illumination. Thus, for a large dielectric surface, 
the current forcing function will have the general form 


Id(V) 





(cos a) 




( 20 ) 


where Iq is the total positive current into a large dielectric surface, A is the area 
of the surface, and all other terms have been defined previously. The above equa- 
tion must satisfy the following condition 

/ 

e®V/X , 


1 If s = +1 ahd V > 0; otherwise leave unchanged 
I if s = -1 and VsO; otherwise leave unchanged . (21) 


Exposed metallic parts of the structure can be located on many different outer 
surfaces of the spacecraft; consequently, the various exposed metallic surfaces, 
which are electrically connected, tan have different sun-lllumlnatlon conditions 
and the current forcing function will have a more complex form. In addition, the 
exposed metallic surfaces are generally small in area. Fewer charged particles 
will be deflected from a small surface at a given potential than a largo surface at 
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the same potential consequently, a correction factor must be applied to small 
surface ar ea s. Thus, the positive eurhent Hewing into the exposed metallic struc- 


ture Is 


-V/T„ / -V/Tj 




t Jpj, (coso^)e 


where Eq. (21) holds for the above equation, is the total exposed metallic 

area, A^ is the exposed area of the l**' metallic surface, m is the total number 
of exposed metallic surfaces, o. is the sun aspect angle for the metallic sur- 
face, and the following holds for the-^mall area correction terms 


(1 + V/Tp) = 


(1 + V/Tp) for V >0 


1 for V < 0 


(1 + 1 V/T I) for v<0 


(1 + 1v/tJ) = 


1 for V > 0 


Equations (20) and (22) are the plasma and photoemission generated current sources 
and constitute the forcing functions of the spacecraft charging equations. 

2.2 Material Properties and ConfipuratloS 

The spacecraft outer surface material '• properties and conflguratlm definition 
are needed to complete both the geometrical and electrical models. Essentially, 
the material properties and configuration definition consist of describing the loca- 
tion of the various outer surface thermal blankets and coatings and their electrical 
properties. The location of the materials is needed In the geometrical model to 
establish the number of constituent planar and curved surfaces of the spacecraft. 
The electrical properties of the materials are needed in the electrical model to 
determine the equivalent circuit element values of the outer surfaces of the 
spelcetfdh* 


The electrical proper^tlea of the outer surface thermal materials that are most 
Impoi^tant H a spacecraft charging aialysU ares 

(1) The relative dielectric constant. 

(2) The vaHatloh of the sUr*faee resistivity with respect to electrical stress 
level. 

(3) The variation of bulk resistivity with respect to electrical Stress level. 

(4) Ratio of surface to bulk leakage currents. 

All of the above properties can be determined experimentally, in fact^ for mean- 
ingful results, the last three parameters should be measured under conditions 
similar to those experienced in the charged particle environment at synchronous 
orbit during a substorm. That is, the measurement results iVill be somewhat 
dependent on the energy levels, and current densities of the charged particles bom- 
barding the dielectric surface of the thermal blanket materials. In practice, how- 
ever, these properties are.measured by bombarding the materials with a mono- 
ergetic electron beam. 

2.3 (k^omelric'al Model. 

The purpose of the geometrical model is to define the spacecraft outer surface 
areas in terms of approximate plariar and Curved surface projections, establish a 
reference coordinate system in the spacecraft, and define the vertices of all of the 
planar and cUrved surfaces in terms of the reference coordinate system. PuHher- 
more, the results of the geometrical model are needed to complete the electrical 
model. That iS, the approximate geometrical surfaces of the spacecraft outer con- 
figuration are used in the computation of the equivalent capacitor and resistor 
element values of the electrical model (each value is related to the surface 
area). The reference coordinate system can be selected anywhere inside the 
vehicle structure and should be chosen such that one or more coordinate axes arc 
parallel to the axes of cymmetry, or parallel to the major planar outer surfaces. 

The reference coordinate system is useful in determining the relative locations 
and orientations of the constituent outer surfaces. In addition, the reference co- 
ordinate system is needed to determine the location of the spacecraft with respect 
to the earth and sun. The surface Vertices are used in the computation of the 
surface normal vectors, and the surface normal vectors together with the location 
of the spacecraft with respect to the sun are used to determine the sUn illumination 
condition of the surface; that IS, whether the surface is illuminated by the sun, 
self-^shadowed, or earth-shadowed. It should foe noted that there can exist outer 
surfaces that are never illuminated by the sun; these areas are designated as 
"permanently** shadowed areas. 


247 


s. 


in order* to deteritilhe the sun*mutnlnation condition of the constituent sur- 
faces all outer surface areas must be expressed in terms of the six major planes 
parallel to the coordinate axes. Hence, the areas of the constituent planar sur- 
Les parallel to a coordinate axis can be easily expressed in terms of the s X 
major planes. However, for planar surfaces not parallel to a coordinate axis 
effective surface areas projected Into the. six major planes must be determine. 

The projected areas in the six major planes are then assumed to have the same 
normal vectors associated With the six major plartOs When determining their s - 
illumination condition; however, the true, unprojected area is used when determin- 
ing the area that is sUl^ect to the charged particle environment. The same proce- 
dure is foUowed for all curved Surfaces, that is. cones, spheres, cylinders. . 

The exposed metallic parts of the spacecraft structure require special 
attention. Since the exposed structural parts can exist on 

planar surface of the spacecraft, the exposed metallic parts Will have differen 
sun-lllumlnation conditions depending on the particular location of the exposed 
part. . The effective projected area of each exposed metallic part is computed in 
Lch of the six major reference planes as outlined above. The effective projected 
area in conjunction with the particular sun/spacecraft surface aspect angle is used 
in the computation of the phOtoemission current; this Is represented by the last 
term of Eq. (22). However, the actual exposed area of a metallic part is subject 
to the plasma envl.onment; thus, the total exposed metallic area is Use^ in the 
computation of the incident particle currents. This is represented by the first two 

t^rms of Eq. (22), 


2. t Eleclficttl Modtl 

The electrical model defines the lumped element, equivalent circuit represen- 
tation of the spacecraft outer surfaces with respect to the electrostatic charging ^ 
phenomenon. The equivalent electrostatic circuit is a network consisting of capaci 
tors and reslstofS whose values are either computed or measured. It wUl be 
assumed that a dielectric surface can he represented aS a simple lumped Capacitor 
and a parallel leakage resistance; however, this Is an approximate representation 
when considering the complex processes that occur when a dielectric surface Is 
bombarded by high energy particles. The capacitor components represent the 
capacitance of the various dielectric surfaces with respect to the spacecraft struc- 
ture The resistor components represent the leakage current from the dielectric 
surfaces to the spacecraft structure. Additional capacitors and resistors are 
needed to represent the surface capacitance and leakage current between adjacent 
surfaces and between illuminated and nortlllUmlrtated sections of a surface. How- 
ever, these surface interaction processes are second-ordfer coupling effects and 
will not be considered in the model. This Is a conservative assumption and does 




not affect the ability of the model to predict the potential differences between a 
surface and the structure or the potential differences between adjacent surfaces. 
Cohsetjuently, the equivalent circuit of the spacecraft with rCspect to the Charging 
phertoitienort has the simplified form shown in Figure 2. 


SURFACE 

elements 


PLASMA 

SOUBC*' 



Figure 2. Spacecraft Equivalent Circuit 


It has been assumed that there are n outer surfaces. The i-th surface has an 
absolute potential of V. volts and each surface, or node, has a. corresponding 
plasma and photoemission generated current source having the general form of 
Eq. (20). The spacecraft structure has an absolute potential of V volts and I is 
the plasma and photoemisslon generated current source into the exposed metallic 
surfaces and is given by Eq. (22). The capacitance, C^, is the intrinsic capaci- 
tance of the spacecraft structure with respect to the plasma. This structural 
capacitance can be approximated by the isolated capacitance of the structure. This 
is a reasonable approximation since the plasma sheath outer boundary, which rep- 
resents the terminus of the strong satellite fields due to spacecraft charging, has 
a depth on the order of tens of meters. 




The follOwlhg set of slmuUeweous spacecraft charglrtg equations caft be written 
foi* the Simplified circuit of Figure 2: 


C, d (Vj 

-V 

(Vi-V„) 

dt 


- V 


-V 

• 

• 

dt- 


- V 




= w 


dV^ 

dt 


n 

S'. 

i-0 


(25) 


EiiuatiDn.(25) in genferal will be nonlinear since the leakage resistances are non- 
linear functions of stress level (V^ - V^) and the plasma and photoemisslon gen- 
erated currents are nonlinear functions of absolute potential. The number of equa- 
tions, n, is a function of both the number of surfaces with different dielectric 
materials and the number of surfaces with different Sun-illUnlination conditions. 


2.5 Solar T^arlh Sp&reiTaft Orbit Model 

The purpose of the solar/earth/spac6draft orbit model Is to determine the 
sun-llluminatlon Condition of a spacecraft surface including both Carth-shadowlng 
and self-shadoWing conditions. The sun-lllumlnation condition of a surface is 
determined by first defining the planar surface and its vertices with respect to the 
spacecraft reference coordinate syst.^m. This Is essentially accomplished in the 
geometrical model. Next, the normal vector of this surface is computed and the 
relative location of the normal vector with respect to the spacecraft reference 
coordinate system Is determined. The relative position of the sun with respect to 
the earth is computed as well as the relative position of the spacecraft with respect 
to the earth, using coordinate transformations, the relative position of the space- 
craft coordinate system with respect to the sun is then determined. Finally, the 
angle between the surface normal and the sun vector, the aspect angle, is com- 
puted. The intensity of sUn-lllumlrtatlort is propOrtlottal to the cosine of the 
aspect angle with full illumination occurring for an aspect angle of 0®. The surface 
is self-shadowed when the absolute value of the aspect angle exceeds 90°, Also, 
the earth-shadowing condition, which occurs when the spacecraft is in the umbra 
of the earth, cart be similarly determined. The solar/earth /spacecraft orbit 
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iiukIpI oonsi-sts of four sepafate parts; the spacecraft ephemeris model, the solar 
model, the coordinate transformation?, and the solar/vehicle/earth geometrical 
model. 

Th^ relative geometry bet^Veen the earth, the sun. the spacecraft, and a con- 
stituent surface IS shown in Figure 3. As indicated in Figure 3, the surface is 
defined by the vertices A. R, and C. To determine the solar/spacecraft surface 
aspect angle, the surface normal vector and the surface-sun vector must be 
computed. The vertex vectors of the surface expressed in vehicle Coordinates 
are 


FA 

FB 

FT 


(26) 


and the sun vector in inertial coordinates at the center of the earth is 


and is computed by the solar model program. The spacecraft vector expressed in 
inertial coordinates is 


OF 


(28) 


and IS computed by the vehicle ephemeris program. 

From Figure 3. it can be seen that the surface normal vector is given by 

N (FA - fT?) w (FC - FR) 

the solar 'Spacecraft surface aspect angle is then given by 
cos o (N • RS),(! Nj IRs! ) . 

and the earth aspect angle is 

COSO (RS - RO) (’Rs! ’rT)') . 
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Ftgurfe 3. Sun-Earth-Spacecrdft Orbital deom^try 


Earth -ecllbsing (shadowing) of a surface is given by the following condition 


3 < 3^ eOltbee Condition 
3 > 3 j non-eclipae condition 


where 


3^ = sin“" (ft^/r) 

= radius of earth 
e 


arid Pj is the earth disc aspect angle subtended the spacecraft. Thus, an earth 
eclipsing condition occurs when the spacecraft is in the umbra of the earth. Solf- 
shadowtng or self-eclipslrig of a surface by the spacecraft itself is given by the 
following condition: 


t 


la I 2r 7 t/ 2 feclipae condltioh 
la I < ff/2 non-eclipso condition 


C^4) 


where a is the solar/spaeecraft surface aspect angle determined from the solar 
projection upon the surface normals. 

In the simplified spacecraft ephenieris model, the location of the spacecraft 
with respect to the inertial coordinate system of the earth is determined. There 
is no need for a precise spacecraft orbit so an abbreviated model is used. The 
orbit is assumed to be circular with a constant radius and a nominal period of 
1440 min. The spacecraft is flowrt in the equatorial plane (inclination = 0^). 

In the solar model, the position of the sim with respect to the earth is derived 
from the American Ephemeris and Nautical Almanac which provides mean position 
in terms of a series expansion of elapsed centuries from a base epoch. Conven- 
tional Coordinate transformations are employed to determine the position of the 
sun with respect to the spacecraft (vehicle) reference coordinate system. 

2.6 Nunicriral Intt'^rUiioti 

The first order, simultaneous, nonlinear spacecraft charging differential 

Eqs. (25) Writh forcing functions represented by Eqs. (20) attd (22) are of such a 

form that standard closed-form methods of solution do not apply. Consequently, 

several "initial value" numerical Integfatlort techniques were utilized to compute 

the time response of the absolute potentials, V. V . The greatest success 

1 1 ^ • 

was achieved with the Runge-Kutta integration process. After some preliminary 
experimentation, it was found that a step size of 0.001 sec produced satisfactory 
results in the time response computation. The step size is the incremental value 
of the independent variable, time, at which the dependent variable value, absolute 
potential, is computed. 

From initial computations of the time response of the set of Eqs. (25) using 
actual circuit values of capacitance and nonlinear resistance and actual plasma 
substorm parameters, it became apparent that the transient and steady-state 
results could not be obtained in a single numeHcal integration execution. First, 
it was found that the steady-state vnlUes of potential are feached after several 
hoUi's. Second, the computet* execution time-to-^solUtion time was enormous 
(typical rill* times were on the ot*der of 20 to 30 min to obtain 1 to 3 min of simu- 
lated time). Consequently, it was decided to characterize the transient behavior 
by computing the transient response up to that point in time at which the transient 
response was well-behaved, that is, either mohotonically decreasing or increasing 
(usually on the order of 1 to 3 min). The steady-state solution was computed 



I I 


seoafately in a fapldly executed progfam. This overall approach was not rlgor- 
oZly accurate alhce plasma conditions can change within minutes and illuminatton 
conditions can change within tens of minutes. However, If worst case plasma and 
Illumination conditions are employed, the solutions, both transient and steady- 
state will represent worst case values and more accurate solutions should not be 


TO obtain the steady-state solution, a more direct method was employed. 
The steady-state condition is characterized by the condition dV^/dt - 0. From 
Elq. (25), this results in the steady-state system of equations 



C: 


(V, - Vq) 

h - (R,(V, - V ) 


1 S i £ n 


and 



^i 


(35) 


(36) 


The solution to this system of equations can be viewed as an optimisation 
problem where Eq. (36). which represents the current balance condition, must be 
minimized while simultaneously satisfying the set of n nonlinear equations. Eqs. 
(35) which can be considered as constraint equations on the current balance con- 
dition. For simplicity, it was decided to use a simple direct enumeration scheme 

to iteratively search the region 


(38) 

V. = V, + rt^v , 0 < I < n 

1 Li 

such that Eq. (36) was minimized while satisfying the rt constraint Eqs. (35). The 
terms V. and V„ are the lower and upper bounds, respectively, of the absolute 
potentials of the surfaces and structure. t>ositlve potential values were incremen- 
ted by aV = 0. 1 volts and negative potehttal Values by AV - 50 volts. 
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Spacecraft charging simulation results for a three-axis stabili/,cd spaciH^raft 
are presented in this section. The spacecraft analyzed, V^lth thermal blankets m 
place, could be adequately modeled geometrically as a ”box-like'’ structure with 
large "flat panel" type solar cell arrays which arc located above and below the 
north and south panels, respectively, of the main spacecraft structure. The 
antenna structures, with thermal blankets, could be modeled as conical structures 
that protrude from the front side of the spacecraft main structure and point towards 
the earth*s surface. The dielectric properties of the thermal blankets and surface- 
coating materials ^^ere measured and the equivalent capacitances and leakage 
resistances of the constituent surfaces were computed. The results are listed in 
Table 1. There were 13 surfaces with either different materials or different 
orientations (with respect to the spacecraft reference coordinate system) that had 
to be considered in the spacecraft charging analysis. The front side had three 
different materials and the north and south panels had two different materials. 

The resistor values listed in Table 1 are based on the bulk resistivity character- 
istic and represent the values Computed at low stress level. The last element in 
the table is the structural capacitance and w'as computed by using some of the 
formulas listed ill Appendix A. 


Table 1. Element Value Summary of Three-Axis Stabilized Spacecraft 
Analyzed 


1 . 

2 . 

3. 

4. 

5. 
n. 

7. 

8. 
9, 

10 . 

lU 

12 . 

13. 

14. 


Element Location 

Backside 
West Panel 
North Panel 
North Panel 
South Panel 
South Panel 
East Panel 
Froht Pahel 
Froht Side 
Front Side 

Solar Array SUh-Side 

Solar Array Dahk-Slde 

Permanently Shadowed 
Sides 

Spacecraft Structure 


Resistor Value (ohms) 

Rj 8. 9 y 10^ 

R 2 1.2 y 10® 

Rg 2. 1 y 10® 

R^ - 1.8 y 10® 

Rr 1.4 y 10® 

r!! 4.0 V 10 *° 

R^ 1.2 V 10 “ 

Ii„ = 2. 5 y 10® 

Rg 9. 5 y lO^'^ 

H,h 2.8 y 10 “ 
Hjj 3.0 V 10^'^ 
Rj 2 1.4 y 10® 

R ” 3.8 10® 


Capacitor Value (pf) 


Cj 0.37 
C 2 0.29 
C 3 - 0.1(5 
c]j 0.1(5 
C 3 0.24 

cj. 0.08 

C 7 0.28 

Cy - 0. 13 
Cj, 0. 009 
C^K 0.024 
Cjj 0.65 
Cj2 4.4 
C,3 - 0.087 

Cj 4 0.000350 
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In genfel’ai, the bulk fealsllvity ia a functtort uf atfCaa level. The bulk feala- 
tivitiea of all the dielectric macCHala were mettsUfed by bombarding .samples of 
the ihateriaLs by high energy electcoas and measuring the through conduction leak- 
age current as a function of electron accelerating beam veltage. ’t M/as Conserva- 
tively assumed that the sur face stress level was approximately equal to the beam 
voltage and a piecewise approximation to the bulk resistivity versus beam voltage 
characteristic was computed. All of the piecewise approximations of the dielectric 
materials had a form similar to the piecewise approximation of Chemglaze palrtt 
shown in Figure 4. To simplify the simulation and to decrease the execution time, 
the piecewise approximation of all of the materials were employed in the analysis. 


/ 



BfiAM VOtTAGi: INkV 

Figure 4. Piece-Wise Approximation of Chemglaze Paint Bulk Resistivity 
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I’ablo 2, Photo^miaaloh, Secondary Emis^^^lurt, arid Omnidirectional Plarima 
Pa^*amdtera 


Parameter 


Typical PangC 


Value Selected 


lV^Tp^-3V 

2 V V 

s 

0 £ f si 
®M 

0 S f <1 
®D 

0 £ f £ I 

Pm 

0 £ f £ 1 

Pd 


0.82 na/cm'^ < Jp 

< 4 na/cm2 

0. 02 na/cm^ < 

< 2 na/cm^ 

2 pa/cm2 < - 

32 pa/cm ^ 


( 6. 0 kV (substorm) 
j 3,0kV(quicl) 
y 20, 0 kV (severe substorm) 


( 12.0 k\ 
{ Vy. 0 k^ 

( 40.0 k^ 


kV (substorm) 
kV (quiet) 

kV (severe substorm) 


2. 0 na/cm2 


0, 6 na/cm2 (substorm) 
0. 02 na/cm2 (quiet) 

0.02 na/cm^ (substorm) 
2. 0 pa/cm2 (qui^t) 


current densities or energies. Consequently, transient solutions were obtained at 
the ortset of the plasma substorm, where particle temperatures (energies) and cur- 
rent densities change suddenly from their quiescent values to their substorm 
values; at the beginning of earth-eclipse, where the photoelectron current forcing 
function is zero; and at the end of earth-eclipse, where the photoelectron current 
forcing function becomes nonzero. The transient response of the structure, the 
solar array, and the surface that exhibited the greatest steady-state potential dif- 
ference Is shown in Figure 5 for the onset of the 6 kV substorm. In this figure, as 
well as the others to be presented, the transient response is shown for a period 
of 70 sec and the steady-state solutions are shoivn on the right side of the figure. 
The transition period fro n the transient to the steady-state solutions is Indicated 
by the dashed lines. (As expected, the surface that consistently exhibited the 


I 
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Flgui-e 5. Transient Response at Onset of a 6 kV Substorm. Fall-Equinox 23:00 LT 

greatest stfeady-state potential difference between the structure and the surface 
itself, was covered with dielectric material that experimentally had the highest 
value of bulk resistivity. ) It was assumed that all potentials were initially at zero 
volts. It can be seen that the surfaces "fall" instantaneously to a few hunured volts 
with small potential differences between the two outer surfaces and the underlying 
structure. This behavior was typical of all of the surfaces of the spacecraft. The 
absolute potentials than "fall" monotonlcally negative until, after a long period of 
time, the final steady-state values are achieved. The transient response at the 
onset of earth-eclipse is shown in Figure 6. It was assumed, as a worst case, 
that the steady-state values of the previous period had been achieved at the start 
of the earth-eclipse and are the initial values used in the numerical integration 
program. It cah be seen that the structure instantaneously "falls" to a negative 
value of about 9 kV. however, the Initial potential differences are maintained but 
decrease mohotonically in the steady-state to small values on the order of a few 
hundred volts. The transient response at the end of the earth-eclipse period is 
shown in Figure 7. Again it Was assumed, as a worst case, that the steady-state 
values of the previous period had been achieved at the end of the earth-eclipse 
pehlod and these values then became the initial values in the numerical integration 
program. The structural potential instantaneously decreases to a negative value of 
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Figuce 6. Transient Response into Eclip'e for a 6 kV Substorm. Fall-Equinox 
23 {30 LT 


a few hundred volts; however, the potential differenees are maintained initially and 
then increase monotonlcally to somewhat larger values. The solar array is 
"clamped" to zero volts (actually less than one volt positive) by photoemission. 

This condition occurs for most surfaces with full sun-illumination intensity. 

Upon comparing the steady-state with the transient solution Values, it becomes 
apparent that the final steady-state values, that is, those values achieved if the 
sun -illumination conditions did not change, represent the worst-case differential 
values. Also, the differential poter^^als can change instantaneously by no more 
than a few hundred volts. But, thr absolute potential of the spacecraft structure 
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Trsuisietit Response Out of Etlipse for a 6 kV Substorm. Fall-Equinox 


ean change almost instantaneously since it has small capacitance. In all Cases, 
the transient response times are controlled by either the potential constants of 
the forcing functions or the time constants of the circuit elements or both. The 
differential potentials of the outer surfaces do not change instantaneously by large 
amounts since the circuit element time constants, ^hich are large in value (the 
product of resistance and capacitance), are dominant. However, the absolute 
potential of the surfaces, Which is the sum of the absolute potential of the structure 
and the differential potential between the surface and the structure, can change 
instantaneously in conjunction With the structure. This is demonstrated in all of 
the transient responses and in particular in Figures 6 and 7. At the beginning and 
at the end of earth-eclipse, the absolute potential of the structure changes by a 
large amount and the absolute potentials of the surfaces change by a similar amount; 
thus the potential differences* do not change IH value initially. 


At the oaset of the eubstorUi, the stfnicture " falls" Slowly negatively since 
there are photoelectton emlSsloha Which tend to cancel the ihfluk of electrons from 
the iJlasma current Source. At the OnSet of earlh-ecUpse there are no photoelec- 
tron currents and the plasma electron Current sources dominant and rapidly 
charge-up the small structural capacitance. At the end of earth-ecllpse, the large 
photoelectron^ current sources again reoccur. The large outflux of electrons from 
the exposed metallic parts produced by the photoemission currents iS instantan- 
eously supplied by the structural capacitance -and consequently, there Is a rapid 
decrease in the negative absolute potently', of the structure. 

As discussed previously, foh three-axis stabilised spaceccaft there is con- 
siderable diurnal as well as seasonal variations In the amount and location of the 
outer surface areas of the spacecraft that are exposed to the sim. TfiUS, the sun 
Uliuninatlon condition of the 13 dielectric sxirfaces of the spacecra^are computed 
throughout the orbital path. From the solar /earth/spaCecraft orbit model It was 
found that the sun -Illumination condition did not change significantly In local time 
increments less than 30 min. Consequently, steady-state solutions were Computed 
at 30 min increments throughout the duration of the substorm. In general, the 
fl'ial steady-state solution will never be achieved at the end of the 30 min period 
slice some of the source potential and network, time constants Involved ace on the 
order of thousands of seconds and the Initial Sim -illumination conditions, on which 
the final steady-state solution IS based, will change significantly every 30 mla. A 
Summary of the spacecraft steady-state values using the average plasma parameter 
values of the 2 January substorm, is listed in Table 3 for 1 Hr Increments through- 
out the duration of the substorm for the fall-equinox period. The hourly Incre- 
mental values are representative of the worst-case potential differences obtained 
when compared against the values computed in the smaller half-hour increments. 
The widely varying values of the absolute potential of the spacecraft structure are 
shown as well as the maximum surface differential potentials. It can be seen that 
during eclipse (lasting roughly from 23:30 to 0:45) the Spacecraft, in steady-state, 
achieves a negative potential Of about 0 kV and a maximum surface potential differ- 
ence of -4. 5 kV was reached towards dawn. 

A similar analysis was conducted for the peak of the winter-solstice period for 
the 6 kV test substorm. As expected, the lowest negative spacecraft structural 
potential was achieved during this period. A summary of steady-state values is 
presented in Table 4. A maximum potential difference of -4. 0 kV was achieved. 

As Can se seen from the summary tables, the spacecraft structural potential 
varies widely reaching a maximum negative value of about 9 kV at eclipse and a 
minimum negative value of 450 kV during the winter-solstice period. In this 
particular design. Ibe structural exposed metal was kept to a minimum and this 




















help^ U> explain the fact that the structure never achlCv^ed ?.erc> potential when 
various exposed parts were illuminated by the siift. Examination of the steady- 
state values of all of the surfaces indicates that during eclipse all of the surfaces 
achieve almost the same absolute potential. This results from the fact that all 
surfaces have the Saine shadow and plasma current source conditions* 


u u)\c:u SIGNS AND DlSCtSSiGN 

This paper has been concerned with the development and application of a 
charging model for three-axis stabilized spacecraft. The objective of the model 
is to determine the differential potentials between the outer surfaces and the struc- 
ture of a spacecraft throughout its geosynchronous orbit when under the influence 
of a geomagnetic sUbstorm. It w'as assumed that the interaction between the plasma 
and the spacecraft can be adequately represented by an equilibrium theory approach. 
That is, the energy distribution jof t!\e constituent plaSma particles Can be expressed 
in terms of an omnidirectional Maxwell-Boltzmarin distribution. The plasma is 
then represented by equivalent voltage dependent current sources and the outer 
Surfaces by simple lumped elements. The resulting first order differential equa- 
tions are integrated and potential distributions determined. Sun-illUmlnation con- 
ditions were determined by a solar/earth/spacecraft orbit model and the intrinsic 
capacitance of the spacecraft with respect tc the plasma sheath is approximated by. 
its isolated capacitance, spacecraft charging simulation results, including both 
the transient and steady-state Solutions, have been presented. 

A .knowledge of the potential distribution of the outer surfaces and structure of 
the spacecraft throughout its orbital path is important from a systems design and 
analysis viewpoint. In general, spacecraft materials that maintain stress levels 
below their dielectric strength level should be selected. If, from the analysis, it 
appears that the dielectric strength of various spacecraft surface materials will be 
exceeded, then, depending on the magnitude and repetition rate of the discharge 
and location of the material, corrective action such as modification or replacement 
of the material may be necessary. Thus, the spacecraft charging simulation re- 
sults can be useful in determining the selection and location of the type of outer 
surface thermal blanket oi^ coating materials to be employed in the design of 
spacecraft. 
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«I«iUUtton results can bo useful in establishing the 

swceoraft charstag simulation reaulla. It was found that when a maximum 

lllumlnaled by the sun, the slnictufal potential 
achleued Us lowest negative value (beeause Of Photoemisslon). At the same time 
the stmfaoe differential potential attained Its maximum value. Conversely when' 
a minimum amount of exposed metal was Illuminated, the spaeeeraft stricture 
echoed tts highest negative value and the surtaoe differential potential attained 
Its rnm^urn value. The above results demonstrate that It I, dLlral to emX 

sinoe the d:reXxxt“:r:,:r r r r*" 

Uhenwm depend on the spaoeeraft eonftguratlon, o'rblt, and olr Xf"'^ 

In the simulation analysis employing a step function, that Is. tlme-lmlepend- 
ent, representation tor the plasma substorm, the worst case dllle.-entl.l potential, 
oceu^ at steady-state and dot during the transient response. ThTIZI^T 

^ widely varying funotlons of time. Because of the large time constants of the 
eoulvalent spaeeeraft circuit, a steady-state response using the actual time 

xrrtTh-“ 

Mle. e thnt the magnitude of the differential potentials as well as the absolute 

iXXt «’» “agnltude of the substorm particle energies 

(temperatures) aild current densities. ^ 

Upon examining the spacecraft charelne results u keah. 
trajectories, their actual energy distributions, and determine the surface 
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pot6ntials by sblving Poisson's potential distribution equation in three dimensions. 
Although the model is based on simplistic assumptions, the simulation results 
Obtained for the structural potential are in relative agreement with the structural 
potentials measured on-board ATS-G, a three-axis stabilized spacecraft. The 
charging model predicts that upon entering eclipse, the spacecraft structure falls 
almost instantaneously to a value of about -3 kV artd leaving eclipse the spacecraft 
structure rises almost instantaneously to a few hundred volts negative. Similar 
transient results into andjbut of eclipse have been observed on ATS-6. In addition, 
the structural potential variations during the post-eclipse period, as predicted by 
the model, correspond, relatively, to measured results on ATS-6. For example, 
examining the steady-state stress levels in the post -eclipse period, it can be seen 
that the structural potential rises to a low negative value after eclipse but falls 
negatively towards dawn and then rises to a low negative value at dawn. These 
results are in relative agreement with data measured on board ATS-6. 
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L IMIiODLiCTlUS 

lA this appendix, Equations for the isolated electrostatic capacitance of several 
types of geometrical structures are given, 

LI Sphere 

The Isolated capacitance of a Sphere is' 


^ISO " 


where is expressed in farads, R is th6 radius of the sphere in meters, and 

is the. permittivity of a vacuum (8.85 x lO'^^ farads/meter). 

K2 Cube 

2 

The capacitance of art isolated cube has been found to be • 

^ISO " ^’^^o ^ ^^2) 

where i is the length of the sides of the cube in meters and is expressed in 

farads. 


l.j (Alindor 

The Capacitance of an isolated cylinder is given by** 


47re_a 


r 2 

a (a^ ^ R^) 
R 


where is expressed in farads, a is one -half the length of the cylinder in 
meters, and R ia the radius of the cylinder in meters. 


This formula was derived bv the authors. 


I <id<* 


I 


i 


The capacllanee of a truhcated cone wUl be approximated by the capacitance 
of a cylinder with a length equal to that of the cone, but the radium of the equivalent 
cylinder Id the average of the radii of the truncated cone. Udlng Eq. (A3), the 
results are 

^"iso " “ 

In 


. 2 *2 
a ^ (a^ ^ R^) 

R 


where H - ^ISO Expressed In farads and the radii expressed in 

meters. 


1.5 Thill Kortan^ailur und t-llliptioal riutos 

The capacitance of thin rectangular and elliptical plates have been derived in 
a previous paper^ and the results are given in graphical form for various values 
of length and width, and semimajor and semimirtor. axes, respectively. 


1.6 Thin i'irruhir IMuir 

4 

The capacitance of a thin circular plate is given by 

=— (Hi. 1) (A5) 

ISO 7T 

where R is the radius of the circular disk in meters and is in picofarads. 
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